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The light pseudoscalar particle, dubbed the axion, borne out of the Peccei-Quinn solution to the 
strong CP problem in QCD remains elusive. One avenue of inferring its existence is through its 
coupling to electromagnetic radiation. So far, laboratory experiments have dedicated all efforts to 
detect the axion in the mass range 10~ 6 < m a < 10~ 3 eV with a photon-axion coupling strength 
g ajJ < 10 -10 GeV -1 , where the limits are derived from astrophysical considerations. In this study, 
we present a novel way of constraining g all by looking at the level of linear polarization in the 
radiation emerging from magnetic white dwarfs (mWDs). We find that photon-axion oscillations 
in WD magnetospheres can enhance the degree of linear polarization. Observing that most mWDs 
show only 5% linear polarization, we derive upper limits on g aT y for different axion masses. 

PACS numbers: 14.80. Va, 97.20.Rp 



I. INTRODUCTION 

Quantum chromodynamics (QCD) has emerged as 
a phenomenologically accurate theory that describes 
strong interactions among the six quark flavors that are 
bound into two families of hadrons, namely mesons and 
baryons. From experiments, we understand that strong 
interactions enjoy C (charge conjugation), P (parity), T 
(time reversal) discrete symmetries of nature. Therefore, 
QCD must also obey such symmetries, both separately 
and any combinations formed thereof [T]. However, CP 
symmetry is broken in QCD due to the presence of the 
following term in the QCD Lagrangian [2] 

£int = ( ^ tr G^Ga^v (1) 

where 9 is a periodic parameter, g is the QCD cou- 
pling constant, G M „ is the color field strength tensor, 
and G^y is its dual. The value of the ^-parameter is 
not set theoretically, but it can be measured from the 
electric dipole moment of the neutron (d n ), for which 
many theoretical estimates exist but we only quote one, 
\\d n \\ ~ 2.7 x l(T 16 #e cm [3]. Here 9 = 9 + arg det m q , 
where m q is the quark mass matrix. The latest exper- 
imental estimate of ||d„|| < 2.9 x 10 _26 e cm [3] con- 
strains \\9\\ < 10 -11 pp. This inexplicably small value 
of 9 gave rise to the strong CP problem. One of the solu- 
tions, also the most favoured, to this problem was envi- 
sioned by Peccei & Quinn [5], whereby the 9 parameter 
is driven precisely to zero under a global chiral symme- 
try, later named £/(1)pq. The pseudo-Nambu-Goldstone 
boson that results upon the spontaneous breakdown of 
this symmetry was dubbed the axion [51 [5] . Not unlike 
the Higgs boson, the axion has proven to be extremely 



difficult to observe as it couples only weakly to ordinary 
matter and radiation. 

Despite several attempts to experimentally observe the 
axion, it remains elusive to this day. Nevertheless, the 
experimental efforts have not gone in vain, but have been 
able to place serious constraints on the coupling strength 
of the axion to photons g ail < 10~ 10 GeV -1 . Stringent 
constraints have been placed on the mass of the axion 
10~ 6 < m a < 1CP 3 eV with the lower limit arising from 
cosmology [7] and the upper limilQ from the neutrino 
flux recorded for SN 1987A, which placed strong limits 
on the cooling flux through other channels namely, right- 
handed neutrinos or axions [TO]- If #077 > 1CV 10 GeV -1 , 
the production of axions through the Primakoff process 
will significantly alter the core He burning timescales of 
post main sequence stars, a possibility excluded by the 
ratio of horizontal branch stars in globular clusters |11) . 
Several reviews on the properties of axions have been 
forthcoming in the past decade, for example see [U fl2| 
113] , to which we point the reader for a more detailed and 
comprehensive exposition. 

Still, there is no denying the fact that none of the lab- 
oratory experiments conducted thus far have been able 
to secure a positive detection of this mysterious particle. 
The detection of very weakly coupled particles demands 
extremely sensitive laboratory experiments. So far, only 
a handful of experiments, namely the Cern Axion Solar 
Telescope (CAST) [HJ, Axion Dark Matter Experiment 
(ADMX) [J3 [T|>], and Rochester-Brookhaven-Fermilab 
(RBF) collaboration [T71 [18] have been able to surpass 
the astrophysically derived limits on <? Q77 in the above 
quoted axion mass range. Yet, the sensitivity envelope 
needs to be pushed even further by a few orders of magi- 
tude to be able to draw any definitive conclusions about 
the existence of the axion. Plans are afoot to modifiy 



1 Due to large uncertainties in the axion mass derived for the DFSZ 
* |rsgill@phas.ubc.ca| model M from SN 1987A observations (0.004 < m a < 0.012 
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the existing experiments and devise new ones to improve 
upon current limits on g ail (see Section IV). 

Unlike the laboratory experiments, the odds are in 
favour for detecting axions in astrophysical systems. This 
optimism stems from the fact that the axion to photon 
conversion probability scales with large magnetic field 
strengths and longer coherence lengths [15] . such that 
■P -> 7 oc g 2 B 2 L 2 , where L is the length over which both 
the photon and axion fields are in phase. Thus, there is 
a very good chance of finding the axion in strongly mag- 
netized compact objects, namely magnetic white dwarfs 
(mWDs) and neutron stars (NSs). The possibility in the 
latter case has been expounded by many (see for exam- 
ple [T9l422j : also see [23] where constraints on <? Q77 are 
derived from the dimming of radiation by photon-axion 
conversion in astrophysical sources), however, the case of 
the mWDs has not been investigated in great detail and 
warrants further study. 



A. Magnetic white dwarfs 

After the discovery of the first mWD by Kemp 24J, 
the number of white dwarfs with magnetic fields rang- 
ing from a few kG to 10 3 MG has grown to about 170. 
The size of this subpopulation is only 3% of the total 
population of known WDs comprising of 5447 object^ 
The main channel for identifying magnetism in WDs is 
through Zccman spectropolarimetry, which not only al- 
lows one to discern the strength of the field but also the 
direction of the field lines, and also through cyclotron 
spectroscopy (see for e.g. [25] for a review on isolated 
and binary mWDs). Nevertheless, reconstruction of the 
field topology has proven to be very difficult, mainly due 
to its highly non-dipolar structure. Over the last decade 
Zeeman tomography of mWDs has enjoyed some suc- 
cess in elucidating the underlying field structure. This 
technique is based on calculating a database of model 
spectra, where different field geometries comprising of 
single/multiple dipole, and higher multipoles, that may 
also be off-centered and misaligned with the rotational 
axis, are considered. Then a least-squares fit using the 
pre-calculated synthetic spectra is performed through a 
highly optimized algorithm on the phase-resolved Zee- 
man spectra to obtain the complex field structures [26] . 
The generality of the models not only allows greater flex- 
ibility but also renders a closer fit to the actual field ge- 
ometry of the source for a given rotational phase. 

The presence of even a small degree of circular polar- 
ization in the spectrum of a WD is a strong indicator of 
the object possessing a magnetic field upwards of 10 6 G 
[2"7] . The degree of circular polarization typically reaches 



2 G.P. McCook and E.M. Sion, web ver- 
sion of the Villanova White Dwarf Catalog, 
:/ /www. astronomy. villanova.edu/WDCatalog/index. html 



up to ~ 5%, and sometimes beyond that in a few se- 
lective objects, near absorption features and also in the 
continuum. Continuum circular polarization stems from 
the magnetic circular dichroism of the atmosphere, where 
the left and right circularly polarized waves propagating 
through a magnetized medium encounter unequal opac- 
ities [3H]. A relatively higher degree of circular polar- 
ization also appears near the red and blue shifted wings 
of the Zeeman split absorption lines {a + and cr_ compo- 
nents, [25]). 

On the other hand, most observations of mWDs indi- 
cate that the linear polarization component never exceeds 
that of the circular one, and the spectrum remains dom- 
inantly circularly polarized until field strengths > 10 8 G 
are reached [28] , In a magneto-active plasma, the plane 
of linear polarization undergoes many Faraday rotations, 
an effect that arises due to the magnetic birefringence 
of the medium, so that on average the degree of linear 
polarization of the emergent radiation is much reduced 

m - 

The very fact that linear polarization is of the order of 
a few percent (~ 5%) in the continuum spectra of most 
mWDs can be exploited to draw meaningful conclusions 
on the extent of axion interaction with photons traversing 
the magnetized plasma of mWDs. We explain how this 
can be implemented in the next section. 



B. Plan of this study 

The purpose of this study is to conduct a survey of 
the m a — g ail parameter space by modelling photon- 
axion oscillations in the magnetosphere of a mWD. To 
this end, we model the field structure of a strongly mag- 
netized WD PG 1015+014, for which high resolution op- 
tical spectropolarimetric observations are available [30] . 
In the same article, the authors also conduct a phase- 
resolved Zeeman tomographic analysis and derive a best- 
fit model of the magnetic field topology. Despite fitting 
the spectrum with a range of field geometries, they were 
only able to pin down the field geometry for a single ro- 
tational phase by fitting it with a superposition of three 
off-centered and non-aligned dipoles of unequal surface 
field strengths (see Table [I] for model parameters). To 
model the effect of photon-axion oscillation in the mag- 
netosphere on the emergent polarization, we propagate 
an unpolarized photon of a given energy from the pho- 
tosphere through the encompassing magnetosphere, that 
has been populated by a diffuse, cold ionized H gas. The 
emergent intensity and polarization is then averaged over 
the whole surface of the star. Finally, we compare the de- 
gree of polarization from our model simulation to what is 
observed in mWDs with field strengths in excess of a few 
10 6 G, for example PG 1015+014, and draw conclusions 
on the strength of the coupling constant for a given axion 
mass. 

In the following Section, we formulate the key equa- 
tions describing the interaction of the axion with pho- 
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TABLE I. Top: Magnetic field geometry adopted from the 
spectropolarimetric analysis by [30] of the the mWD PG 
1015+014. The model comprises of three off-centered and 
non-aligned dipoles D\,D2,D% with unequal surface field 
strengths B s , polar (9b) and azimuthal (4>b) angles of the 
magnetic field axes. The center positions of the dipoles rela- 
tive to the center of the star are given by (a x , a y ,a z ). Bottom: 
Collection of parameters used in the study. 



Model Parameters 



D-l D 2 D 3 


B s (MG) 


-40 92 -38 


0b{°) 


44 63 63 


4>b(°) 


339 276 134 


a x (R+) 


0.04 -0.012 0.27 


a y (R*) 


0.35 -0.136 0.080 




0.33 -0.28 0.21 


R* 


7 x 10 s cm 


0, 


23° 


Y e 


1 


Bq 


4.413 x 10 13 G 


T 


10 4 K 


9* 


10 8 cm s" 2 


Po 


10~ 10 g cm" 3 


Poo 


10~ 20 g cm' 3 



tons, geometry of the aggregate magnetic field, and struc- 
ture of the plasma permeating the magnetosphere. The 
lack of understanding of the density profile of the magne- 
tospheric plasma introduces some level of inaccuracy in 
any treatment of mWDs. We take the simplest approach 
and describe the plasma density by the barometric law 
for an isothermal atmosphere. In Section III we present 
the main results of the study along with a comparison to 
some of the results obtained from lab experiments (see 
Fig. [5]) . A discussion of the results is provided in Section 
IV. 



II. MODEL EQUATIONS 



Z 




X 

FIG. 1. This illustrates the coordinate system used to obtain 
the aggregate magnetic field. Here Eo represents the coordi- 
nate system centered on the star with the z-axis aligned with 
the rotational axis. Ej represents the coordinate system in 
which the different dipole magnetic field equations are writ- 
ten. This system is misaligned with Eo by a polar angle 9b,% 
and an azimuthal angle (j>B,i, and then it is displaced from 
the center by the vector a^. For clarity we have chosen the 
vector a.i to lie along the z'-axis (color online). 



coupling strength, E is the polarization vector of the pho- 
ton field, and B is the external magnetic field. Quan- 
tum corrections to the classical electromagnetic field, due 
to the constant creation and annihilation of electron- 
positron pairs in vacuum, are given in the last term of 
([2]) by the Euler-Heisenberg effective Lagrangian, causing 
the vacuum to be birefringent. 



The interaction of the axion field with an external elec- 
tromagnetic field is given by the following Lagrangian 
density |20| . in natural units where h = c = 1, 

£ = - iF^F^ + \{d^a - m 2 a a 2 ) - ig^F^F^a 

2' 



'Mm 1 



(F^F^f + 1 (f^p^ 



(2) 



A. Off-centered non-aligned three dipole model 

We start by writing down the aggregate magnetic field 
Bo in the coordinate system Eo where the z-axis coin- 
cides with the rotational axis of the star 



where the first term describes the external electromag- 
netic field, with F^v as the antisymmetric electromag- 
netic field strength tensor and F^ v = ^e^ l/pcr F pcr as its 
dual. The second term is simply the Klein-Gordon equa- 
tion for the axion field a where m a represents its mass. 
The next term is the interaction Lagrangian density, 
which upon simplification, using the given definitions, 
yields £i n t = <? a77 aE • B, where g a -yy is the photon-axion 



(3) 



The three off-centered dipole fields B^ are first rotated 
before they are added together by operating on each of 
them with a rotation matrix R T = R^, i R^ )il where the 
superscript T indicates the transpose, and the rotation 
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matrices are given as follows 



Rz.i 



Ry,i — 



cos(f> B ,t sm<f> B ,i 
-smc/) Bti cos (f>B. i 
1 

cos 6>b,i — sin 0b,., 

1 
sin 9b i cos 9 B .i 



(4) 
(5) 



Here the polar and azimuthal angles 9s,i and 4>B,ii re- 
spectively, are defined with respect to the axis of rotation. 
In cartesian coordinates, the dipole fields are expressed 
as 



B' 



Bs.iR 



±(3^x' + 3y^y' + {32f -rf}z') (6) 



' 2 i 



2r' 



where the fields are shifted from the coordinate center, 
such that r' i = r — a,j. In the above equation, B s i is the 
surface field strength of the i th dipole field component, 
~ 7 x 10 8 cm is the radius of the WD, and r\ is the 
magnitude of the radial vector in the coordinate system 
Ej. The profile of the aggregate field B as a function of 
distance is shown in Fig. [2j 



B. Fully ionized pure H atmosphere 



The presence of a magnetic field necessarily introduces 
anisotropy in the plasma dielectric tensor s^. In the case 
of a nonuniform field, none of the dielectric tensor com- 
ponents vanish, as compared to the homogeneous case. 
Below we write all the dielectric components, which one 
can easily derive from Maxwell's equations, for complete- 
ness. 



£ll £12 £13 
£21 £22 £23 
£31 £32 £33 



(7) 



~|2 d2 
J c,s a Qx 



1 



•,2 



l-£c,e(l + £c,p)£oa 



£12 



£13 



£21 



£22 



£23 



(1-^e) 

- 2 

~(i-^L) 

- 2 

"(i-fL) 



(l-^ e )(l-^ P ) 



" e ;{iB 0z + Ld c e B Qx B 0y ) 



(iB 0y — 6j c ^Bq x B 0z ) 



(iBqz — &c,eBo x Bo y ) 



p,e 
;2 



1 - u>l a (l + u>l p )B { 



(iB Qx + Lo c ^B Qy B Qz ) 



(i-aL) 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 



£31 



' — 3r^(iB 0y + Lo c ^ e B 0x B 0z ) 



£32 ~ ~ ' ; V \ {iBQx — £>c,eBo y B() z ) 



£22 ~ 1 - w 



l-^ c 2 e (l+u; c 2 p )i3 2 , 
(l-<tf e )(!-<&» ) 



(14) 
(15) 
(16) 



In the above set of equations, w c s = q s BQ/ujm s c is the 
normalized cyclotron frequency for species s = (e,p), 
where e and p signify electrons and protons; w PjS = 
\/ 'ATTn s /m s w 2 is the normalized plasma frequency, where 
Tip = n e = Y e p/m p are the electron and proton number 
densities, Ye is the electron fraction, and p is the pro- 
ton mass density of the plasma; the normalized magnetic 
field components are defined as So,i=x,j/,z = Bq^/Bq. 



1. Vacuum corrections 

Due to the polarizability of the vacuum in strong mag- 
netic fields, the plasma dielectric tensor g^, and the in- 
verse permeability tensor p~ x are modified |3lJ|32], such 
that e —e = e + Ae and = = I + A11" 1 , 

-p+v - -p -v h=p+v = = =v 

where 



Apr 1 = (a v — 1)1 + m„BoB 

a v = 1 — 2S V q v = 7S V m v — —4S V S v = 



(17) 
(18) 



-Bo 

457T \ Bq 



and Bq = 4.413 x 10 13 G is the quantum critical field for 
which the separation in energy between Landau levels of 
the electron exceeds its rest mass. 



2. Plasma density profile 

That many mWDs are surrounded by hot coronae has 
been suggested by many to explain the polarized flux of 
those WDs that show comparable degree of linear and 
circular polarization [29l |33l [34] . The thermal electrons 
in the hot tenuous plasma with temperature T ~ 10 6 ~ 8 K 
radiate at the cyclotron frequency that falls in the optical 
wavelength for field strengths of B ~ 10 s G. This radia- 
tion appears to be polarized both linearly and circularly, 
depending on the orientation of the line of sight to the 
magnetic field, and traverses the corona without any ap- 
sorption. Furthermore, slightly polarized radiation em- 
anating from the photosphere, with very low degree of 
linear polarization due to Faraday rotation, gets added 
to that generated in the corona, as a result increasing 
the amount of flux that is polarized linearly. Several hot 
isolated WDs, with effective temperatures in excess of 
~ 25, 000 K, emitting X-rays were detected by ROSAT 
[35] . however all cases were linked to subphotospheric 



5 



thermal emission [Mj . Although the non-detection of any 
coronal emission may indicate the absence of a hot ten- 
uous corona, it is not at all unreasonable to suggest the 
presence of a tenuous cold plasma of fully ionized H. In 
this study, we envisage that the mWDs are encompassed 
by cold isothermal electron-proton coronae with the fol- 
lowing barometric density profile, 



p(r) = p exp 



(19) 



where po is the density near the surface of the star, p x 
is the density that remains far away from the star as 
the strength of the magnetic field becomes significantly 
weaker than that at the surface, and H p = 2kBT/m p g* ~ 
1.65 x 10 4 cm is the density scale height with an ef- 
fective temperature T ~ 10 4 K and surface gravity 
log(?*(cm/s ) = 8. There is no clear agreement on the 
surface plasma density with 10~ n < po < 10~ 6 g cm~ 3 . 
Here, we assume that the plasma is sufficiently tenuous 



with po = 10 g cm 3 and p Q 



= 1(T 20 g cm- 3 . 



C. Axion-photon mode evolution in an 
inhomogeneous magnetized plasma 

We are interested in knowing the evolution of the axion 
field and the polarization vector as the radiation prop- 
agates out from the surface of the star, traversing the 
region with an inhomogeneous plasma density and mag- 
netic field. Here we follow the discussion given in [20]. 12 1] . 
and derive the photon field mode evolution from the EM 
wave equation 



V x (fr 1 • V x E) = f • E 



(20) 



Next, we assume the ansatz E = Eexp(ifcz) where the 
wave is propagating along the rotational axis of the star, 
which in this case is also the line of sight direction, and 
the wavenumber k — uj/c. Plugging this ansatz into the 
wave equation, and ignoring second order derivatives, we 
find 



d 
dz 



E. r 

E„ 



Xn Xl2 
X21 X22 



E u 



-*y j y A.^i a.^^ y y ^y 

where the matrix elements are given below 

T 2 



(21) 



Xii = V 



k'iu - T 4 - 1 - 



t,t. 



TiT i 



x I k 2 e 2 i 



T 2 



^ 1 r; 2~ 
L 3 



T 4 } 



(22) 



C12 = V 



T 2 

TsTf 



^£12 - T 2 - ^1 - 

k 2 e 22 - T 6 - ^{k 2 e 12 - T 2 } 
13 



TiT 1 



(23) 



X21 = 1 - 



Tf 



I k 2 e 21 - T 2 - ^{fc 2 en - T 4 } 



(24) 



X22 = 1 - 



T 2 
TaT^ 



T _i 





/ Ti 

x (k 2 e 22 -T 6 -^-{k 2 e 12 - 


T 2 }J 


(25) 


Ti = 


d - - - 
— {m v B x B y ) + i2km v B x B y 




(26) 


T 2 = 


ik^(m v B x B y ) - k 2 m v B x B y 




(27) 


T 3 = 


— -r(a v + m v B 2 ) + i2k(a v + 
dz " 


m v B 2 ) 


(28) 


T 4 = 


— ifc — (a„ + m v B 2 ) + k 2 (a v - 


-m v B 2 y ) 


(29) 


T 5 = 


~~dz^ av + m "^) + i^k(a v + 


m v B 2 x ) 


(30) 


T 6 = 


—ik — (a v + m v B 2 ) + k 2 (a v - 


- m v B 2 ) 


(31) 



1. Line of sight geometry 

The Zeeman tomography analysis of mWD PG 
1015+014 indicates that the line of sight (LOS) is in- 
clined at an angle 0k = 23° to the rotational axis of the 
star. Following [37], we modify the matrix Eq.plj) to 
obtain the mode evolution of the photon-axion system in 
a coordinate system oriented along the LOS (see Fig. |2j. 
Again, we assume the ansatz a oc exp (ik' s — iujt) 



. d 

ds 



where A n 



a \ 




( K-k' 


Ami' 


&My' 


E x > 






i X 'n - k' 


ix'12 


Ey, J 




V Am v > 


1X21 


ix'22 - k 



a 

E x i 

Ey' 

(32) 



to 2 /2u;, A M x> = -g ail B x /2, A My > 
—g ail B y /2. Notice that Eq.plj) applies to a system for 
which the LOS vector coincides with the rotational axis 
of the star. For a different LOS vector, such as shown in 
Fig. [2] we perform a rotation of the plasma dielectric ten- 
sor around the y-axis by an angle 8k, — Ry {6k)eR y (9k) 
where R y is given in Eq.([5|. 

The total degree of polarization can be found by in- 
tegrating Eq.(32) from a given point on the surface 



outwards to a distance beyond which the amplitude of 
photon-axion oscillations and plasma effects become neg- 
ligible, and then by averaging the Stokes parameters [55] 
over the whole observable hemisphere. 



/ = 



\E x i 
\E X > 



\Ey\ 
\E,A 



(33) 
(34) 
(35) 

-i(E x ,E;,-E v ,E*,) (36) 
The mode amplitudes are in general complex, and in 



U 
V 



E X ,E* 



Ey'E*, 



6 




ur 

lo 6 
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10" 
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FIG. 2. This figure illustrates the coordinate system used to 
obtain the photon-axion mode evolution along a given LOS 
(top), and the decline of the magnetic field strength with dis- 
tance s from the surface (bottom). Here the LOS vector is 
represented by s that is tilted at an angle 9k to the rotation 
axis, and is in parallel to z' which is not to be confused with 
z' in Fig. [I] Several different points on the star's surface with 
spherical coordinates (R*,9u,(j)u) are chosen and then aver- 
aged to determine the final polarization of the photon leaving 
the magnetosphere (color online). 



the above set of equations * gives the complex conju- 
gate. We sample the Stokes vector from different points 
on the surface, with coordinates (R+,9r,4>r), that are 
spread around the LOS vector k with A(f> R = 30° and 
A0 R = 10°, where 10° < R < 80° and 0° < (j) R < 330°. 
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FIG. 3. Polarization evolution of an unpolarized photon along 
a given LOS as it starts at the photosphere and propagates 
through the magnetosphere. In this case, Ey = 3 eV, m a = 



10- b eV, g a 



lCT a GeV . The magnetic field geometry 



assumed is that of mWD PG 1015+014 (color online). 



Because the sampling in the azimuthal angle is sparse for 
larger polar angles, we take a weighted average, as shown 
below for one of the Stokes parameters, to determine the 
average degree of polarization of the whole hemisphere 



(I) 



(37) 



III. RESULTS 

In the following, we look at how an unpolarized photon 
emitted from the photosphere of a mWD gets polarized 
as it traverses through the magnetosphere. Photon-axion 
interaction and the intervening plasma make the medium 
birefringent, consequently, altering the state of polariza- 
tion of the unpolarized photon. We obtain the degree of 
polarization from the averaged Stokes parameters 



Pl 
Pc 



(I) 



(I) 



(38) 
(39) 



where Pj, and Pc represent linear and circular polariza- 
tion. 
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In Fig. [3j we present the evolution of the Stokes vec- 
tor with distance s from the surface of the star for the 
case of radiation with £L = 3 eV, and axion parame- 
ters m a = W- 5 eV,g a77 = 1(T 9 GeV -1 . The oscilla- 
tions in the solution arise due to the mixing of the ax- 
ion and photon eigenstates, an effect analogous to neu- 
trino oscillations due to the MSW effect [39l|40]. How- 
ever, notice that the interaction is non-resonant because 
a 50% drop in intensity would be observed if the ax- 
ion and photon modes were to achieve maximal mixing 
and undergo level crossing. Eventually, as the photons 
travel farther away from the surface, the decline in the 
magnetic field strength reduces the probability of con- 
version, hence the diminishing of intensity variation. We 
hnd that the change in polarization is primarily brought 
about by the axion interaction with the photon. In the 
event this interaction is made negligible, no significant 
polarization or change in intensity of the emergent ra- 
diation is found. The origin of circular polarization in 
mWDs, as alluded to earlier, is understood in terms of 
the difference in opacities for the two modes of radiation, 
making the plasma dichroic, in the presence of a mag- 
netic field. Linear polarization, on the other hand, was 
explained by the cyclotron radiation that emanates from 
the tenuous corona composed of an ionized plasma. In 
this study, since the treatment of radiative transfer ef- 
fects is very simplistic only an upper limit can be placed 
on how strongly the axion couples to photons, as shown 
in the next section. 
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FIG. 4. The final state of polarization of radiation after 
traversal from the WD's magnetosphere for different m a and 
g a -yy- Here (7) is the average Stokes intensity, and Pl and Pc 
are the degrees of linear and circular polarizations. These re- 
sults apply to the case of mWD PG 1015+014 (color online). 



A. Constraints on g all 

Axion production in the mWD magnetosphere can en- 
hance the degree of linear polarization of the observed 
optical radiation. The goal here is to not determine the 
precise value of the photon-axion coupling strength but 
only constrain it from above. To this end, we look at 
the amount of linear polarization that is produced for a 
given m a and g 077 - The underlying assumption here is 
that all of the observed linear polarization is generated 
due to photon-axion interaction, and not by the plasma, 
which effectively yields the absolute upper limit on c/ a77 . 
In Fig. |4j we plot the emergent intensity and state of po- 
larization for different axion masses and for photons in 
the UV - optical waveband with energies between 2 — 5 
eV. The m a and <? a77 in Fig. [3] were chosen specifically 
so that P L > 0.05 for all photon energies. 

In Fig. [5] we use the same parameters to draw an ex- 
clusion region in the m a — g Q77 parameter space, along 
with regions excluded by lab experiments and astrophys- 
ical considerations. The shaded region in red excludes 
all m a - g ajl values for the case of mWD PG 1015+014, 
that is for a typical surface field strength B ~ 10 8 G 
and degree of linear polarization Pj, ~ 5%. We find that 
for the range of masses that are of relevance, in partic- 
ular, to the axion models, the constraints on <7 Q77 from 
this study are superseded by that from horizontal branch 



(HB) stars. Still, the limiting linear polarization criterion 
used in this study is able to probe smaller <7 a77 values in 
comparison to works that only look at radiation dimming 
(for e.g. see [23]). 

The constraints can be further improved by looking at 
mWDs with higher magnetic field strengths. The highest 
field strength that has ever been discovered in a mWD is 
B ~ 1000 MG in two such objects namely, PG 1031+234 
and SDSS J234605+385337 gj. Both objects show lin- 
ear polarization as low as ~ 1% for some rotational 
phases [JH H3J. Based on these two facts and assum- 
ing that the magnetic field geometry of these two mWDs 
is at least as complex as that found in PG 1015+014, we 
produce two exclusion regions shown in Fig. [5] with col- 
ors blue and green. The former corresponds to a surface 
field strength B = 1000 MG with the same level of lin- 
ear polarization as before, and the latter studies the case 
with Pl ~ 1%. For these two cases, we have only looked 
at m a < 10 -5 eV since higher mass values don't con- 
strain <7q 77 better than limits derived from HB stars and 
CAST (Phase-I). On the other hand, we have extended 
our treatment to smaller particle masses with < 10~ 6 
eV where should be interpreted as the mass of any 
light pseudoscalar boson that is characteristically very 
much similar to the axion but isn't a CDM particle. 

It is worth mentioning that the change in <? 077 is not 
linear with the change in magnetic field strength, as ev- 
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FIG. 5. Exclusion plot in the m a — g all parameter space. The 
red region corresponds to the case of mWD PG 1015+014 
with B ~ 10 s G and limiting linear polarization Pl ~ 5%. 
The blue and green regions correspond to the case with 
B ~ 10 9 G, but with P L ~ 5% and P L ~ 1% respectively. 
The mass of the axion is constrained to 10 -6 < m a < 10" 2 
eV from cosmology and SN 1987A measurments. The photon- 
axion coupling constant is capped from above with <? a77 < 
10 -10 GeV -1 by the number of horizontal branch stars in 
globular clusters. KSVZ ggj [47] and DFSZ E] are two 
different theoretical models that predict how g a -,-, scales with 
m a . Other exclusion regions are from lab experiments by 
the CAST experiment [14], ADMX group [HI [16] and by 
Rochester-Brookhaven-Fermilab collaboration |17l 118] (color 
online). 



ident from the comparison between the red and blue re- 
gions in Fig. [5] For higher field strengths one observes 
higher degree of polarization of the emerging radiation. 
Naively, one would expect the plane of polarization to 
rotate by an amount that is 0(g 2 71 B 2 l 2 ), which is valid 
strictly in the absence of plasma when the photon and ax- 
ion are treated as massless particles [501 SI]- Therefore, 
for a fixed degree of polarization an increase in B should 
also decrease g ail by the same factor, when I, the length 
over which the magnetic field remains homogeneous, is 
kept constant. However, as shown by [45] in the case 
of NSs, an increase in magnetic field strength also in- 
creases the level of polarization by effectively shifting the 
polarization-limiting radius, the distance beyond which 
the two polarization modes couple and which depends 
weakly on the magnetic field strength R p i oc B 2 / 5 , farther 
away from the star. The farther the polarization-limiting 
radius, the more coherently the polarization states from 
different LOSs add, yielding a higher degree of polariza- 
tion. 



IV. DISCUSSION 

This study looks at how the production of axions in 
mWD magnetospheres can alter the state of polariza- 
tion of the observed radiation. We find that unpolar- 
ized photons of photospheric origin become linearly po- 
larized upon their traversal through the inhomogeneous 
magnetic field of a mWD. We have modeled the magneto- 
spheric plasma, as fully ionized pure H with a barometric 
profile. Since the majority of mWDs are strongly circu- 
larly polarized and only show a relatively small degree of 
linear polarization, at most 5%, we have used this obser- 
vation to constrain the coupling strength 5077 of axions 
to photons. We find that for the case where the plasma 
component only contributes negligibly to the state of po- 
larization, the coupling strength .9,277 increases with the 
mass of the axion m a - The level of linear and circular 
polarization observed in mWDs is sensitive to the prop- 
erties of the magnetospheric plasma. The limits on gayy 
can be improved by modelling all the radiative transfer 
effects in WD atmospheres and fitting the model spectra 
to real observations. 

Magnetic fields stronger than that of mWDs exist in 
NSs. Going back to the argument of how astrophysi- 
cal objects, compared to laboratory experiments, ben- 
efit from longer coherence lengths (see Section I), in 
comparing mWDs with NSs, one finds that the latter 
are ~ 10 4 times more efficient in converting photons to 
axions and vice-versa. A number of studies have ex- 
pounded on the subject of propagation of polarized ra- 
diation through the NS magnetosphere, where they have 
considered IR/Optical radiation [38], and thermal X-rays 
[32 HH] produced at the surface of the NS. Unfortunately, 
no X-ray polarimetry observations have been conducted 
partly due to the very low flux in X-rays from these 
objects, and also because none of the high energy tele- 
scopes are equipped with a polarimeter. X-ray polarime- 
try has been neglected for the last 30 years but it is hoped 
that some of the future space missions [50 , for exam- 
ple the Gravity and Extreme Magnetism Small Explorer 
(GEMS) [51] . will fill this void in X-ray astronomy. In 
any case, as discussed by [T9H21"] . NSs are excellent lab- 
oratories for the detection of any light, weakly coupled 
pseudoscalar particle. 

A. Outlook 

The ADMX[^] project, that employs a microwave cav- 
ity to search for cold dark matter axions, will begin its 
phase II of testing in the year 2012. With the new up- 
grades the ADMX project will be able to exclude g ajl 
up to the DFSZ line in the same mass range as before. 
Although outside of the range of axion masses probed in 



3 http:/ /www. phys.washington.edu/groups/admx/experiment. html 



9 



this study, the modified CAST experiment has been able 
to exclude axions with g ail > 2.2 x 1CT 10 GeV -1 for 
m a < 0.4 eV, becoming the first experiment to ever cross 
the KSVZ line [SI [S3]. The currently running CAST 
experiment in its phase-II will be able to exclude axions 
with m a < 1.15 eV with unprecedented sensitiviy in this 
mass range. An improved version of the light shining 
through wall (LSW) experiment ([33], also see [S3] for 
a recent review on such experiments) using Fabry-Perot 
optical cavities to resonantly enhance photon-axion con- 
version has been proposed |55H57| . The projected limit 
in sensitivity to <? a77 > 2.0 x 10~ n GeV -1 typically for 
axion masses m a < 10~ 4 eV achieved using 12 Teva- 
tron superconducting dipoles appears quite promising. 
Further improvements in experiment design and opti- 
mization techniques yielding increased sensitivity to even 
smaller coupling strengths have also been suggested by 
many workers in the field, for example the use of the 
dipole magnets, each providing a field strength of 5 T, 
from the Hadron Electron Ring Accelerator (HERA) at 
DESY in Hamburg in a 20+20 configuration can poten- 
tially exclude g ail > 10" 11 GeV" 1 for m a < 10" 4 eV 
[581 159"] . Another proposed line of investigation to search 
for axion like particles (ALPs) is the use of resonant mi- 
crowave cavities which are much similar in design to the 
optical LSW experiments discussed above [501452] . This 



method has already been employed to search for hidden 
sector photons [63 and can prove to be a powerful tool 
in the case of axions. 

Finally, the simplistic model assumed for the mWD at- 
mosphere only yields an absolute upper bound on g a -yy- 
A much tighter constraint can be obtained by adopting a 
more realistic atmospheric model and solving the equa- 
tions of radiative transfer with the photon-axion oscilla- 
tions included. Such an analysis is outside the scope of 
this study, but it is hoped that the novel method dis- 
cussed in this work will prove to be extremely useful in 
better constraining the properties of any ALP. 
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